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Abstract 

The wavefunction of two fermions, repulsively interacting in the presence of 
a Fermi sea, is evaluated in detail. We consider large but finite systems in order to 
obtain an unambiguous picture of the two-particle correlations. As recently pointed 
out by Anderson, in d < 2 dimensions the particles may be correlated even when 
situated on the Fermi surface. The "partial exclusion principle" for two particles 
with opposite spin on the same Fermi point is discussed, and related to results from 
the T-matrix approximation. Particles on different Fermi points are shown to be 
uncorrelated in d > 1. Using the results for the two-particle correlations we find 
that the orthogonality effect induced by adding an extra particle to a (tentative) 
two-dimensional Fermi liquid is finite. 
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1. INTRODUCTION 

Several years after Anderson |]J conjectured the failure of Fermi liquid theory in two- 
dimensional interacting electron systems even at weak coupling, the issue is still rather 
controversial. Many body perturbation theory says that Fermi liquid theory breaks down 
only in ID, while it seems valid in any higher dimension, at least at weak coupling. 0] 
However, the wavefunction for two interacting particles in the presence of a Fermi sea 
exhibits a peculiar and at first sight alarming feature, namely a finite phase shift for 
two particles with opposite spin sitting on the same point of the Fermi surface. jl| A 
finite phase shift signals the presence of correlations in the two-particle wave function, 
which seems at odds with the existence of independent and stable quasi particles. This 
effect is present at arbitrarily weak coupling in 2D. Anderson [[J indeed argued that due 
to the finite phase shift the insertion of an extra particle in a 2D Fermi system causes 
an orthogonality catastrophe, making the quasi particle weight vanish, and leading to 
Luttinger instead of Fermi liquid behavior. 

Within conventional many body theory the phase shift is blurred in the thermody- 
namic limit by coarse-graining momentum space, and does not seem to have any drastic 
consequences. However, the perturbative many body formalism may be inadequate when 
singular correlations in momentum space appear. Taking the thermodynamic limit before 
solving the interacting problem is dangerous in this case. To obtain an unambiguous 
picture of the correlations associated with the finite phase shift and its possible conse- 
quences, it is therefore worthwhile to go back to Schrodinger's equation, and analyse the 
structure of wave functions in large finite systems. 

In this work we present a detailed analysis of the correlations between two locally 
interacting particles in the presence of a Fermi sea, in one, two and three dimensions, 
extending earlier studies by Anderson Jl| and by Stamp. || In Sec. 2 we provide some 
basic definitions and concepts useful for a clear discussion of large finite systems. In 
Sec. 3 we solve the two-particle Schrodinger equation, and discuss in detail how the 
wave functions are affected by the interaction. A careful derivation of Anderson's "partial 
exclusion principle" |I[ for two particles on common Fermi points will be given, and the 
controversial relation between "phase shift" and "phase angle" {|, [j| will be clarified. 
Particles on different Fermi points are shown to be uncorrelated in d > 1. In Sec. 4 we 
will show that the orthogonality effect induced by adding an extra particle on the Fermi 
surface of a (tentative) Fermi liquid is finite in any dimension above one. Hence, in higher 
dimensions a breakdown of Fermi liquid theory, if any, must be more subtle than in ID. 
Finally, in Sec. 5, we will conclude with a few remarks on the possibility of hitherto 
undetected non-perturbative phenomena at weak coupling. 
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2. LARGE FINITE SYSTEMS 

Our aim is to take the large system limit only after having understood the effects of 
interactions. It is therefore necessary to define all quantities appearing in the course of 
the calculation for finite systems, and introduce certain distinctions which are usually not 
made in the infinite volume limit. 

For definiteness we consider a one-band Hubbard model on a d-dimensional simple 
cubic lattice with lattice constant one and periodic boundary conditions. The Hilbert 
space of states may be spanned by antisymmetrized products of local single particle states 
cjjO), where cj a creates a fermion with spin projection o on site j and |0) is the vacuum. 
Alternatively one may construct a basis from states with sharp momentum |ker) = a k(T |0), 
where a kcr = V^ 1 ^ 2 J2j e% ^ c \en an d V = L d is the number of lattice sites. We assume L to 
be even. The momenta k are taken from the Brillouin zone B — {k = (ki, .., kd) : k v = 
(— L/2 + 1, — L/2 + 2, . . . , L/2)^ 1 }, which forms a lattice with V sites and lattice constant 
2n/L. The Hamiltonian is 

H = H + H I = Y,t»<lc Jff + (2-1) 

ij> j 

where tij is a (translation invariant) hopping matrix, U > a (repulsive) coupling con- 
stant, and rij CT = c^Cjo-. The kinetic part can also be written in diagonal form as 

#o = E e k4 CT a k(7 (2.2) 

k,cr 

where ek is the Fourier transform of iy. 

The iV-particle ground state of the non-interacting system (U = 0) is given by 

\K) = n n °lio> (2.3) 

where the non-interacting N-particle Fermi sea T N contains all momenta in B with < 
fj, N , and the chemical potential fj, N is such that T N contains A^/2 momenta (N = N-^ + N^, 
assume — i.e. N is even). In the following we will frequently drop the index N. 
Let T = B\J r be the complement of the Fermi sea in B. Vectors k^ denote momenta 
in T on the border to JF, while k^ denote momenta in f on the border to JF (see Fig. 
1). The set of all k^ and kj form the "inner" and "outer" Fermi surface dT~ and <9jF + , 
respectively. In the limit L — > oo both sets define the same manifold dJ 7 , the Fermi 
surface. We assume N < L d and e k to be such that the Fermi surface is smooth and 
convex without nesting. 

In the following, when considering very large finite systems or, more precisely, sequences 
of larger and larger systems, it will be important to distinguish various levels of "equality" 
of momenta: 
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i) "microscopic" equality, k'=k, if both are precisely on the same site of the k- lattice, 

ii) " mesoscopic" equality, k' = s k, if both may be separated by a finite number of steps on 
the k-lattice, 

iii) "macroscopic" equality, k' = c k, if both may be separated by an infinite number of 
discrete steps (as L — > oo), which is smaller than O(L), however, such that the distance 
between k' and k shrinks to zero for L — > oo. 

Properties described for k ,n ^k or k' = s k will be understood to hold for "almost" all such 
cases (zero measure exceptions allowed); k'"=k is what is usually implied by writing "k' = 
k " in many-body theory, when performing calculations directly in the thermodynamic 
limit, where "k" refers actually to (infinitely) many states, and it is supposed that they 
need not be distinguished any more, maintaining only their density V/ (27r) d in momentum 
space as the only information. On the other hand, the Pauli exclusion principle acts 
only in the case of "microscopic" equality, but is fortunately easy to build in exactly, 
and is all one needs of k-space fine structure in the non-interacting system. In general, 
in an interacting system it is not a priori clear whether the "internal" structure of a 
"point" in the continuum of momenta in the infinite volume limit is really irrelevant. A 
simple (though admittedly unphysical) example for an interaction where it is relevant 
would be a strict exclusion principle for particles with opposite spin on the same (in the 
microscopic sense) point in k-space. The Hubbard or other short range interactions are 
of course smooth in momentum space, but singularities in k-space might be generated 
non-perturbatively. 

3. TWO PARTICLE WAVE FUNCTION 

In this section the wave function for two particles with opposite spin in the presence of 
a Fermi sea will be evaluated. As in the Cooper problem, the Fermi sea will be assumed 
to be inert, i.e. its role is merely to block momentum space. Much of the calculation 
in (A) and (B) follows the analysis of the related problem of two particles on an empty 
lattice by Fabrizio, Parola and Tosatti.|J 

A) Schrodinger equation: 

The wave function for two particles with total momentum P (conserved) is a linear 
combination 

|*> = E' I P /2 + q T, P/2 - q I) (3.1) 
q 

where the prime restricts the momenta P/2 ± q to T . The amplitudes L(q) obey the 
Schrodinger equation 

(E - £°(q))L(q) = -^E' W) = : C (3.2) 

q' 
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where E°(q) := ep/2+ q + ep/2-q- For U > there are two classes of solutions, a trivial 
class characterized by C = 0, and a non-trivial one with C^O, respectively. 

In the former case one has eigenvalues E = E° where E° is a non-interacting eigenvalue, 
and L(q) ^ only for q such that E°(q) — E°. In addition, the amplitudes are restricted 
by the condition X^L(q) = 0. For each d°-fold degenerate E°, there are d° — 1 such 
solutions, where dP/2 are spin-triplet and d /2 — 1 spin-singlet (if d > 2). Usually d° is 
at least two, due to the symmetry E°(q) = E°(—q), an exception being E°(0). 

In the latter class, one can solve for L(q), and obtains 

L(ci) = ° (3.3) 

& — CP/2+q — e P/2-q 

while the eigenvalues E are determined by 

JJ = JdT! ( E - e P/2 +q - ep/2-q)- 1 =: K L (P, E) (3.4) 

Note that Kl(P, E) is a real function, which has simple poles at the non-interacting two- 
particle levels E°(q). The normalization 1 = (^|^) = X) q ' (-^(q)! 2 °f the wave function 
fixes C as 

C- 2 = 2' (E- ep/2+q - ep/2-q)- 2 (3.5) 
q 

For fixed total momentum P, the non- interacting two-particle levels E°(q) = e P / 2+q + 
ep/2-q can be ordered in an ascending sequence Eq, E®, E^- In terms of {E®}, the 
eigenvalue equation reads 

1 1 M d° 

U = T^ E^E° (3 - 6) 

where dP a is the degeneracy of the non- interacting level E®. The right hand side has simple 
poles in E®. Hence it is obvious that the solutions of (3.6) also form an ascending sequence 
E , E u E M , where E° a < E a < E° a+1 . 

Let |$kk') denote the non-interacting eigenstate obtained by forming the symmetric 
(spin-singlet) linear combination of |k| k'J,) and |k'| k|), and possibly other states with 
the same energy E° = £k + ek'- For each |$kk') there is a corresponding exact eigenstate 
l^kk') of H, related to |$kk') by continuity as U 0. In the following we will analyse 
the energy shift and the modification of these wave functions by the interaction. In 
particular we will calculate the overlap of interacting and non-interacting wave functions 
as a convenient and easy-to-understand measure for the wave function change, alternative 
to the "phase shift", which will also be discussed. 

B) State l^kk) and "partial exclusion principle": 
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Let us now analyse the state |^kk) i n dimensions d = 1,2,3. We will determine the 
energy shift SE = E(U) — E°, the overlap Sk := (^kkl^kk) m the large volume limit, 
and, if this overlap is smaller than one, the "range" of the interacting wave function in 
k-space. Setting P = 2k and extracting the term with q = 0, the eigenvalue equation 
(3.4) becomes 

where SE := E — 2e k and AE°(q) := tk+ q + e k - q — 2e k . Note that, for small q, 

d 

AE°(q) ~ ^2 a u ql , q u = (2tx / L)n v where n v integer (3.8) 

i.e. the smallest AE°(q) are of order 0(L~ 2 ). 2D turns out to be a critical dimension 
here, because the distance from E°(0) = 2e k to the next non-interacting levels -E°(q) is 
of order L~ 2 in any dimension, while the potential energy of |$kk) (as a trial state) is 
U/L d . The overlap is given by the amplitude L(0) = C/5E, where the normalization 
constant C can be written as 

°~ 2 = W? + 5 ISE-AEO^)]* (3 ' 9) 
For E above 2ek but below the next non-interacting level i?°(q), one has 

K L (2k, E) -> (9(1) ind = 3 (3.10c) 

The exception for k = c kp~ in ID is due to the complete blocking of states close to |$kk) 
by exclusion from the Fermi sea T . In d > 1 there is no such complete blocking for 
k = c kp due to degrees of freedom parallel to the Fermi surface. Partial blocking makes 
Kl(21l,E) finite in 2D for k = s kp~, while it diverges logarithmically otherwise. Inserting 
the asymptotic behavior of K L (2k, E) into the eigenvalue equation (3.7), one obtains the 
energy shifts 

iE fO(L->) brkSkf md=1 



0{L- 2 ) fork^k 



tE^\° { "2 , f0rk I; kF ; in^=2 (3.116, 
\C(l/L 2 logL) fork^ s k+ V ; 
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5E -> 0(L- 3 ) ind = 3 (3.11c) 

Rewriting (3.9) as (S k )- 2 = {5E) 2 /C 2 = 1 + (o"£) 2 £ q y[o"£ - AE°(q)}~ 2 , one thus 
obtains the overlap 

} ~l md=l (3.12a) 

for k^k+ 




f r < 1 for k= s ki" , , 

S k -> < n ™ T md = 2 3.126 

\ 1 -0(l/(logL) 2 ) fork^k+ 1 ; 

S^^1-0(L- 2 ) ind = 3 (3.12c) 

For k'^kp in ID and k= s kp~ in 2D, the behavior of Ki(2k,E) implies energy shifts SE 
of the order of the level spacing, L~ 2 , and thus a finite reduction of the overlap Sk, for 
any non-zero interaction U. In ID, the asymptotic overlap r does not depend on U, as 
long as U > 0, while 2D r is U-dependent and goes continously to one for U — > 0. For 
k^'kp in 2D, 5E turns out to be of order 0(1/L 2 logL), i.e. too small for transferring a 
finite amplitude to other states besides |$kk)- For the exceptional case k"=kp in ID, the 
energy shift 8E is of order L -1 , but here the next allowed levels are separated by a gap. 
The reader may compare with the corresponding results for two particles on an empty 
lattice in Ref. 0. 

The overlap reduction Sk < 1 implies that two particles with opposite spin cannot fully 
occupy the same k-state, a phenomenon which Anderson [J] refers to as "partial exclusion 
principle" . To clearly see this effect it was important to take the limit L — > oo only after 
having calculated the overlap for finite systems at finite U. The wave functions l^kk) are 
very short-ranged in k-space: The amplitudes L(q) in |\l/kk) are of order L~ 2 as soon as 
q differs macroscopically from zero, i.e. for L — » oo the wave function |^kk) is confined 
to an infinitesimally small region in momentum space, and is therefore macroscopically 
indistinguishable from the non-interacting state |$kk)- In contrast to the case of Pauli 
exclusion, the state |$kk) is not (even partially) expelled from the Hilbert space of states. 
For L — > oo, there is not even a partial transfer of amplitude to higher energies: Summing 
the squared probability amplitudes l-L^O)! 2 of states with total momentum P = 2k and 
energies E a in an infinitesimal interval around 2ek, the total occupation probability one 
is recovered. 

C) States |^kk') = 

We now analyse the states |^kk') f° r generic momenta k and k'. The energy shift 
5E = E — ek — £k' can be determined by splitting the eigenvalue equation (3.4) as 

v = z4 + L '%' =■ 4e + k ^ + k ' B > < 3 - 13 > 
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and studying the asymptotic behavior of Kl(\h + k', E) for large L. Here d° is the degen- 
eracy of the non-interacting level E° = tk + e^, and -E' (q) := £k+ q + £k'- q ; the double 
prime indicates E°(q) ^ E° in addition to exclusion of k + q and k' — q from the Fermi 
sea T . The overlap S'kk' = (^kk'l^kk') is given by \^dPC/SE, where C is obtained from 

C ~ 2 = + [E-l^W (3 ' 14) 

Recall that |$kk') is a symmetric combination of d° degenerate states (with amplitude 
1/y/cP for each). For L — > oo, the right hand side of (3.14) is always dominated by 
levels in an infinitesimal interval around E. The qualitative behavior of the sum in (3.14) 
follows from the mean spacing AE° of levels around E, which is related to the density of 
two-particle states 

p(E) = lim lim L- d Y!^ E ~ ^P/2+q - ep/2- q ) (3.15) 

q 

via 

AEO = d°(E)/L d p(E) (3.16) 

Here 5 v (x) is a broadened delta-function of width rj, and d°(E) is the level degeneracy as 
determined by symmetry (accidental degeneracies possible for certain dispersion relations 
have zero measure). We will now discuss results for the overlap S'kk' m various distinct 
cases. 

If k or k' (or both) are macroscopically distant from the Fermi surface, and in addition 
k'm^k, a j eve j S p ac ing of order L~ d around E°, and is finite for E between E° 

and the next non-interacting level. Hence 5E is of order L~ d , as the level spacing, which 
implies that finite amplitude is transferred to other levels, i.e. Skk' < 1 in any dimension 
in this case. 

If k n =k' macroscopically distant from the Fermi surface, one finds 



>kk' 



' r < 1 in d = 1 

1 ind = 2 (3.17) 
. 1 in d = 3 



as is easily understood by extending the corresponding results for k = c k' in (3.12). 

Let us now consider the important case where both momenta lie macrospocially on the 
Fermi surface, i.e. d(k, dT), d(k', dJ 7 ) = c 0, where d(.,.) denotes the euclidean distance 
between points or sets in k-space. In this case, the overlap obeys 

ind=l (3.18a) 
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for W m - W anH d((k+k')/2,9^ fi ., 

tor k - k and finite ind = 2 (3.186) 

else 

ind = 3 (3.18c) 

for large systems. In deriving these results, the three cases k^k', k"^ — k' and k"^ ± k' 
must be treated separately. 

For k"^ ± k' (possible only in d > 1, for k and k' on dJ 7 ), the two-particle density of 
states p(k + k', ek + ek') vanishes, i.e. the levels next to ek + are typically at infinite 
distance on scale L~ d . Hence no amplitude is transferred to other levels, and thus S'kk' - ► 1 
for L — > oo. 

For k"= — k' (Cooper channel) the density of states is finite, but Kl(Il + W, E) diverges 
logarithmically for large L (and E detached from non- interacting levels), i.e. the Schro- 
dinger equation forces SE down to order l/L d logL, implying S'kk' = 1 — ^(l/(l°g-^) 2 ) m 
any dimension. 

For k n =k' the density of states is divergent in ID, zero in 3D, and has a rather subtle 
behavior in 2D. Let us discuss only the most difficult (and important) 2D case in detail. 
The general qualitative behavior can be understood by assuming a quadratic dispersion 
relation = k 2 /2 for simplicity. Setting k F = 1, the density of two-particle (or two- hole) 
states in 2D is then given by M 



p{E) 



for uj < Uq 

1 /Air for ujn < uj < uj- 

■ -1 0,/P r ^ ^ 3.19 

sin —i=L — for cj_ < uj < uj + 

1/4-7T for uj + < uj 



Here uj is the energy relative to 2ep, i.e. uj = E — 1, and the various regimes are separated 
by uj + = uj + (P) =P + P 2 /2, uj_ = uj_(P) = -P + P 2 /2 and uj = uj (P) = (P/2) 2 - 1, 
respectively, where P = |P|. See Fig. 2 for an illustration of the various regimes in the 
(P, a;)-plane. Note that here we make use only of the part where uj > 0, corresponding 
to two particles, not holes. This density of states has a simple geometric interpretation: 
For quadratic dispersion, equi-energy manifolds for two particles with total momentum 
P are spheres (in 2D circles) with center P/2 in k-space; in 2D, the density of all states 
(irrespective of whether particles are in T or JF) is l/4n independent of P and uj. The 
density of two particle states in T in 2D is thus simply x/An, where x e [0, 1] is the 
fraction of diameters crossing the equi-energy circle with both ends in T (see Fig. 3 
for an illustration). The same holds analogously for two holes. The singular behavior 
of p(E) in the limit P — > 2k p, E — > 2tp is simply due to the fact that x may assume 
any value between zero and one, however close to the Fermi surface the particles may 
be. Thus it is clear that generically two particles in k and k' with k' = c k find other 
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levels within an energetic distance of order L~ 2 , the only exception being the cases where 
x = 0, corresponding to c£((k + k')/2, dJ-)/d(lL, k') — > 0. It remains to see under which 
conditions ifi,(k + k', i£) is finite. In almost all cases there is an infinite number of levels 
above and below E° = + e^i since we have required only macroscopic equality of k 
and k'. Hence, for E = E° + SE detached from non-interacting levels, we may check the 
finiteness of the sum Kp(k + k.',E) for L — > oo from the corresponding principal value 
integral KeK(k + k', E), where 

K(P, E) := lim lim K L (P, E + irj) (3.20) 

7]— >0 L—>oo 

Note that it doesn't matter whether we insert Kp or Kp on the right hand side. For 
P > 2kp, K(P, E) is just the particle-particle bubble known in many-body perturbation 
theory (for P < 2k p, however, K differs from the bubble since the two-hole contribution 
is absent in K). In the regimes which are of interest here, i.e. P > 2kp and E — 2ep > ujq 
not large, KeK(P, E) can be taken from earlier results for the particle-particle bubble in 
2D 0, i.e. 

f-flog 4(^o)(^) for u < lj < u_ 

ReK(P,E) = I 477 h«+vW-"X w — " )1 2 (3.21) 
[ log - L p2-^ for uj^ < u < u)+ 

where u c is an ultraviolet cutoff. Obviously Rei^(P, E) does not have a unique limit for 
|P| — >• 2kp, E —> 2tp. Generically Rei^(k + k', + e^) is finite in the limit k, k' — > dJ 7 , 
k' — > k, being divergent only (logarithmically) if the ratio d(k, k') /<i((k + k')/2, dJ 7 ) goes 
to zero in the limiting process (note that d((k + k')/2,dJ 7 ) = P/2 — kp and d(k, k') = 
2(u-uj ) 1 / 2 ). 

Hence, generically the overlap Skk' is reduced for k and k' macroscopically on the same 
Fermi point in 2D, exceptions being the rare cases where the ratio d((k + k.')/2,dJ r )/ 
d(k, k') is either zero or infinite. The geometry of the generic and the two exceptional 
cases is shown in Fig. 4. Viewed as a limiting process where k, k' — > dJ 7 , the phase 
space for this wave function modification in the forward scattering channel vanishes with 
the same power as the one for Cooper scattering. Judging from Fig. 2 the effect in the 
forward scattering channel looks weaker because it corresponds to points in a quadratically 
narrowing region in the (P,uj) plain, for u — > 0, while the Cooper processes take place in 
the only linearly narrowing region u > lj + (P). However, the figure shows only a section in 
the d+1 dimensional (P, to) space, where the Cooper processes take place in a cone around 
the cu-axis, while the 2 /cp-pro cesses live in a quadratically narrowing fissure encircling the 
k>axis at a fixed distance 2k p. Hence, in both cases the phase space vanishes quadratically 
in the low energy limit in 2D. 
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d(k, dF) + d(k', dTT= 


d{k,dF) + d{k',dF) in ^O 


d 


k m =k' 


k^k' 


k m =k' 


k^k' 


1 


r (a) 


1 


r 


r 


2 


r (6) 


1 


1 


r 


3 


1 


1 


1 


r 



Table 1: Values of the overlap Skk' : r denotes an overlap < 1; (a) the special case 
k = k ,n =k+ has Sw = 1, (b) S kk , < 1 generically for d((k + k')/2, &F)/d(k, k!) = finite 
7^ and also for k = c k'*=?kp and k= s k' ==kp (including k n =k' = c kp). 



The results for S^k' in d = 1,2,3 dimensions are summarized in table 1. Analogous 
results hold for two holes instead of two particles. 

D) Phase shift versus phase angle: 

The phase shift Xa for two interacting particles in a finite system is defined by [|T], [5| 

Xa = -nSE a /AE° a (3.22) 

where 5E a = E a — E°, A£° = E° a+1 - E° a . Recall that a = 0, 1, 2, . . . , M labels all 
the different non-interacting two-particle energies E® a of two particles with fixed total 
momentum P in a sequence of monotonously increasing energies, and the interacting 
eigenvalues obey E® < E a < E® +1 for a = 0, . . . , M — 1. The phase shift is a measure for 
the modification of the non-interacting wave function |<3> Q ) by interactions. It is finite if 
(& a \^ a ) < 1 and zero if ($q,|\I/ q ) = 1. The term "phase shift" derives from an expression 
of the form (3.22) for the phase shifts in a partial wave decomposition in scattering 
theory. In d > 1, Xa is a wildly fluctuating function of a, which requires a proper 
average over many levels in order to obtain a well defined limiting function for L — > oo. 
The phase angle 4>{E) is defined by 

T(E) = \T(E)\exp[i(j)(E)} (3.23) 

where T(E) is the 2-particle scattering vertex for an infinite system, which is given by 

T(E) = U -— - (3.24) 

and K(E) is obtained from Kl(E), (3.4), via the limiting procedure K{E) = lim^o 
lim^oo Kl(E + in); the dependence on the total momentum has not been written here. 
Note that our K(E) is slightly different from the particle-particle bubble in perturbation 
theory, since in Kl(E) two- hole contributions are absent. 
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Phase shift and phase angle are in general different quantities, even for large L, except 
in d = 1. In M their equivalence has been shown with the tacit assumption that Xa tends 
to a continuous function x{E) as L increases, which is however not generally the case. 
Only in ID x a has a unique limit for L —>■ oo and — > E, and one can indeed show that 

X (E) := lim Xa = <P{E) (3.25) 

and the overlap of interacting and non-interacting wave functions is related to the phase 
shift by the simple identity 

- ^ (3.26) 

in this case. A derivation for these relations is given in Appendix A. 

Phase shifts and phase angles for two particles with fixed total momentum and variable 
energy are shown in Fig. 5 for a ID system and in Fig. 6 for a 2D system. The phase 
shifts have been calculated for a Hubbard model with next neighbor hopping t = 1 and 
interaction U = 5, on a large but finite lattice. In ID the relation (3.25) is seen to be 
verified, while in 2D the phase shifts fluctuate around the phase angle. In 2D one may 
define a function x(-E') representing the mean phase shift obtained by averaging many 
points in a small energy interval. For the systems studied here it turned out that this 
mean phase shift generically differs from the phase angle (at small energies it is larger), 
but x{E) behaves qualitatively similar to <f)(E), i.e. one has 

X (E) = a(E)4>(E) (3.27) 

where a(E) is a smooth function of order one. In particular, a finite phase angle implies 
a finite mean phase shift and vice versa. Furthermore, following the steps in 0, it is easy 
to see that there is a general identity relating the phase angle to the average energy shift, 
namely 

(f)(E) = -n5E/AE° (3.28) 



where AE° is the mean level spacing. 

A controversy in previous studies arose on whether the phase shift for two particles on 
the same point of the Fermi surface in 2D is finite or not, the problem being that the limit 
P — > 2kp, E — > 2ej? is not unique. It was noticed that the asymptotic phase angle is finite 
if the limit is taken in a particular way. 0, |7], [8| However, as we have pointed out above, 
what looks so special a limit in the (P, E) plane reflects actually the generic behavior of 
two particles with momenta k' and k in the limit k' — ► k — ► dT . Anderson's [|I| finite 
phase shift, which has originally been calculated in the special case k = k', is generically 
present for k^k', too. 
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The detailed behavior of the phase angle 0(ek + ek') with P = k + k' for two particles 
in k and k' in the limit k' — > k — > dT is illustrated in Fig. 7. 

E) Antibound states: 

To complete the presentation of the two-particle problem in a lattice model (with an 
upper bound in energy), we now briefly discuss properties of the "antibound state" [[!]] on 
top of the two-particle spectrum. 

Since E a < E® +1 for a < M — 1, all levels but the highest (for a given total momentum 
P) are shifted only by a tiny amount of order L~ d or even less in some cases. However 
SE M = E M — E M turns out to be finite if U exceeds a critical value U c , which depends 
on density and dimensionality. In d < 2, (3.4) implies that U c = 0, since the density of 
two-particle levels is finite (in 2D) or divergent (in d < 2) at the upper band edge. For 
small U, 5 Em is exponentially small in 2D, while it is of order U 2 in ID. This split state 
has been called the "antibound state". [0] It is in a sense the mirror image of the bound 
state in the Cooper problem with attractive U. 

The antibound state \^m) is actually a bound state in that it has a finite expectation 
value for double occupancy of sites in real space even for L — > oo, i.e. 

(Vm\Hi\*m) ~ 0(1) for U>U C (3.29) 

This is of course energetically highly unfavorable for repulsive (positive) U, which is why 
SEm is positive of order one. Since the trace of Hj in the subspace of two-particle states 
with fixed total momentum is always U, independent of the basis, equation (3.29) implies 
that 

M-l 

(^alH^a) <U for U>U C (3.30) 

a=0 

Hence the splitting of an antibound state means that the other states in the continuum 
have an overall reduced expectation value for double occupancy, and \^m) just pays the 
bill for all of them. In the Hilbert space of two-particle states spanned by {l^a) : a — 
0, . . . , M — 1} double occupancy has been partially projected out. In 2D, this projection 
is a weak coupling effect (present for any U > 0), in 3D not. 

Ignoring states which are separated by a gap from the low-energy part of the spectrum, 
one may say that the presence of a down-spin, say, reduces the dimensionality of the space 
of available states for up-spin particles, as in the presence of statistical interactions [|l(J 
between opposite spins. The significance of the antibound states in the two-particle 
system for the many-particle system is however not yet clear. Anderson [J] suggested 
that the splitting of the antibound states implies that the Hubbard model in d < 2 might 
have the same low energy behavior as the tJ-model, where doubly occupied sites are 
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projected out completely. A full projection of double occupancy can be implemented by 
gauge fields, leading to singular effective interactions which have been argued to invalidate 
Fermi liquid theory for the tJ-model (at least at finite temperatures). Thus, if the above 
(controversial!) arguments were valid, Fermi liquid theory would break down in the 2D 
Hubbard model for any U > U c = 0. 



4. ORTHOGONALITY CATASTROPHE 

Long ago Anderson |12| pointed out that a local scatterer in a many fermion system 
changes the wave function so drastically that its overlap with the wave function without 
scatterer is zero in the infinite volume limit, and related this "orthogonality catastophe" 
to the observed singular response of electron systems to a sudden appearance of local 
scatterers, such as in the X-ray problem. More recently he proposed to extend this line 
of reasoning to the insertion of a quasi particle in an interacting Fermi system. 

To understand the argument, it is useful to recall the case of a local scatterer first. A 
system of non-interacting (spinless) fermions in the presence of a local potential on site 
is governed by the Hamiltonian 



^ = E e k^k + f/L" d ^4 +q a k (4.1) 

k k,q 

The ground state of if is a Slater-determinant constructed with single particle wave func- 
tions of the form |\l/ k ) = I] q L k (q)a k+q |0) where L k (q) = C k /(i? k — e k+q ), the eigenvalue 
E k is the solution next to e k of U^ 1 = L~ d J2q(E ~ e k+ q ) _1 , and C k is fixed by normal- 
ization. Note that SE^ = — e k and C k are both typically of order L~ d , corresponding 
to the spacing of the non-interacting levels (except for k^O). The overlap of the Fermi 
seas with and without scattering potential, respectively, is given by [[12] 



S := = det (L kp ) (4.2) 

k^eJ 17 

where L kp := L k (p — k). Note that L kp decays rapidly as a function of e p — e k , but 
the momentum transfer p — k may be large. A sufficient condition for the orthogonality 
catastrophe S — > when L —>■ oo is that the sum 

*■■= E E I^pI 2 (4-3) 

diverges. Indeed s behaves roughly as L~ 2d J^keJ 7 Sp6^( e k — e p)~ 2 > which is logarithmically 
infrared divergent for L — > oo in any dimension. Note that contributions come from any 
q = p — k (not only small ones) across the Fermi surface. To obtain the orthogonality 
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catastrophe it is important that all (or at least a finite fraction of) the single particle 
states ajjO) with k e dT are modified by the scatterer, i.e. ($k|^k) < 1 or, what is 
the same, the phase shift Xk : = SE^/AE^ must be finite all over the Fermi surface (AE^ 
being the distance to the level following ek). 

Recently Anderson [Q] suggested to infer the breakdown of Fermi liquid theory in two- 
dimensional interacting electron systems from an orthogonality catastrophe caused by 
insertion of an extra particle. To this end he considers the overlap 

ZL=\(alK\*t +1 )\ 2 (4-4) 

for "k = kp". Here ^/^ is the exact ground state with iV particles while ^P^ 1 is an 
exact eigenstate of the N + 1 particle system with "one quasi particle added", i.e. the 
state evolving adiabatically from the non-interacting state aj^ | ) , as the interaction is 
switched on. In analogy to the local scatterer problem it is argued that the finite phase 
shift inflicted by the extra particle on the other particles on the Fermi surface will lead 
to an orthogonality catastrophe, Z^ F — > 0, and consequently the elimination of the quasi 
particle peak in the spectral function. 

This argument presents various difficulties which we will now discuss. 

(i) The macroscopic spectral function p(k, £) is given in terms of exact eigenstates by 

p(k,0 = hn l lim E\(*Kn\*l\K)\\(t-Q (4.5) 

' n 

for positive energies £ > 0; the sum runs over all (N+l)-particle eigenstates with mo- 
mentum k and spin a (relative to |^^)), £ n is the excitation energy, and 5 n a broadened 
5-function, e.g. S n (x) = r]/ir(x 2 + rf). Each "point" (k, £) in p(k, £) involves actually an 
infinite number of eigenstates of the interacting system. Hence the vanishing of the over- 
lap with a single eigenstate in the large system limit does not necessarily affect p(k, £). 
In particular, it is easy to see that Z k — > for generic k =kp in any dimension, even in a 
Fermi liquid. By definition, in a Fermi liquid the spectral function obeys the asymptotic 
behavior 

p(k,£)^Z k 5r k (£-£ k ) (4.6) 

for £ — > 0, k — > dJ 7 , where E^ is the quasi particle energy, Z k a finite renormalization 
constant, and r k the width of the quasi particle peak, which must vanish more rapidly 
than the quasi particle energy when approaching the Fermi surface. In a Fermi liquid, for 
k™kp the width of the quasi particle peak is zero on scale one, but generically infinite 
on the scale set by the level spacing. Hence, Z k is zero in this case. To have a chance to 
get a finite Z k , one must set at least k = kp , which would however leave the numerical 
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value of Z'^ F completely arbitrary, if finite (e.g., for k a hundred steps on the k- lattice 
away from dJ 7 , is much smaller than for k = c kp. A reasonable unique definition of 
a possibly finite ZL requires the choice k = c kp in (4.4). We are not aware of a general 
identity relating Z^ p to the renormalization factor Z^ F in (4.6). A priori Z k might still 
vanish even if Z^ F is finite (but not vice versa, of course). 

(ii) In contrast to the local scatterer problem, the extra particle inserted here has 
a complicated dynamics, and the overlap (4.4) cannot be calculated exactly. Making 
an estimate in analogy to the case of a local scatterer added to a Fermi gas amounts 
to making two (independent) approximations, which may miss important physics: The 
system in the absence of the extra particle is treated as non-interacting, i.e. the ground 
state becomes simply a product of two non-interacting Fermi seas for up and down 
spins, i.e. \^q) ~ \&o)- Only interactions between the extra particle (with spin down, 
say) and particles with opposite spin (i.e. up) are kept. Still the calculation of l^k^j 1 ) 
poses a many-body problem, due to effective interactions between up-spins mediated 
by the extra down-spin: an up-spin may scatter the down-spin to a new state, which 
changes its relation to other up-spins. [fL3[l For a local scatterer (without internal degrees 
of freedom) this problem does not occur, because the scatterer remains always in the 
same state. Neglecting these induced correlations, too, one may estimate the overlap 
{(i^pa^o l^kpa) m analogy to the problem of a local scatterer by approximating 

(4 F ,<\^) « det (L kp ) (4.7) 

where L^ p := Lk(p — k), and the amplitudes Lic(q) are extracted from the interacting 
two-particle wave functions evolving from |$kk F ), i-e. 

|*kk F > = E'^(q)|k + q T k F - q |) (4.8) 
q 

solved for fixed k' = c kp and variable k e T in presence of an inert Fermi sea of down-spins, 
but no other up-spins. Here only the down-spin momentum kp — q is blocked by exclusion 
from T . Note that only the up-spins reduce the overlap, since the down spin Fermi sea 
is treated as inert. Since the up-spins are not blocked by a pre-existent Fermi sea, most 
l^kkp) will be shifted from the corresponding non- interacting states |$kk F )- However, 
most modifications inside T in k-space cancel out when the Fermi sea is filled, and only 
shifts leading out of J 7 , as measured by the sum s = Ske.? 7 SpejH-^kpl 2 , are relevant. 
These latter shifts are of the same order of magnitude as those considered in Sec. 3, 
where |\I/kk') has been analyzed for k and k' on the surface of two already pre-existent 
Fermi seas. 



16 



iii) In Sec. 3 we have seen that two particles on the Fermi surface may be scattered 
within T only if they are on the same point of the surface. Analogously, an extra down- 
particle added to the system in kp in the presence of an inert Fermi sea down-spins may 
scatter up-particles out of T into T only in kp itself, but not on other points of the 
surface. Even worse, if kp = c kp, in most cases the extra particle is not able to scatter 
any up-spin out of JF (this is slightly different from the situation in Sec. 3, where both 
particles were outside T from the start, and thus could always shift parallel to the Fermi 
surface in d > 1). Comparing this state of affairs with the response to a local scatterer, 
where a wave function change over the whole Fermi surface led to a logarithmically (only) 
divergent signal of orthogonality, it is obvious that here we find no signal at all. 

In summary, a straightforward adaption of the local scatterer calculation to the prob- 
lem of inserting a dynamical particle into an interacting many-body system does not signal 
an orthogonality catastrophe in 2D. Stamp, || too, concluded that considering finite sys- 
tems within a two-particle scattering approximation does not yield any evidence for an 
orthogonality catastrophe. Clearly, approximating the ground state by a Fermi gas may 
give qualitatively correct results only if the exact ground state is a Fermi liquid. Hence, 
as in perturbation theory, we have only checked consistency of quasi particle behavior as 
a hypothesis. We have to recognize that the above check of orthogonality is insufficient, 
if the exact ground state is neither a Fermi liquid nor a state obtained by resumming di- 
vergencies showing up in perturbation theory (as in the one- dimensional Luttinger liquid 

MM)- 

Note that the phase shift as calculated in Sec. 3 does not signal the orthogonality 
catastrophe which is known to occur in a ID interacting Fermi system upon adding an 
extra particle near one of the two Fermi points! A particle inserted in k = c kp cannot 
kick out any states near kp itself in ID (and for k = s kp it can affect only a finite 
number). Hence interactions with particles near the same Fermi point do not produce 
an orthogonality catastrophe in ID, whether the phase shift is finite or not. On the 
other hand, interactions with particles on the opposite Fermi point do affect infinitely 
many states, and second order perturbation theory does indeed indicate an orthogonality 
catastrophe in this case. However, the phase shift calculated in Sec. 3 turned out to be 
zero in this case, vanishing logarithmically in the large system limit. This is an artefact 
of our treating the Fermi surface as inert, not allowing for particle-hole excitations when 
calculating two-particle correlations. In a diagrammatic language, treating the Fermi 
sea as inert means summing only ladder diagrams, which is equivalent to introducing a 
renormalized coupling whose flow is calculated from the particle-particle channel only, 
and therefore seems to scale to zero logarithmically (for positive bare coupling). In ID 
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there are however other contributions, involving particle-hole excitations, which make the 
/3-function vanish identically, [111, |T^] i.e. the renormalized coupling and the exact phase 
shift in the many-body system remain finite, and an orthogonality catastrophe does occur. 



5. CONCLUSION 

Two fermions in the presence of a Fermi sea can have interaction induced correlations 
even if both particles are situated on the Fermi surface: In addition to the well-known 
Cooper pair correlations for attractive interactions, in low dimensional systems correla- 
tions appear even for purely repulsive interactions, namely i) if both particles sit on the 
same point of the Fermi surface (in d < 2), and ii) if they sit on opposite points (in d — 1). 
In the former case a solution of the Schrodinger equation for a two-particle wave function 
in presence of an inert Fermi sea reveals these correlations, while in the latter case one 
must allow for particle-hole excitations to obtain the correct result for a many-body sys- 
tem. Equivalently, in the former case a properly interpreted T-matrix calculation yields 
the effect, while in the latter a complete one- loop renormalization group calculation of 
the two-particle vertex is required. 

In the full many-body problem we must actually distinguish correlations between bare 
particles and correlations between low energy excitations (i.e. quasi particles in a Fermi 
liquid). Correlations among bare particles are of course always present in an interacting 
theory, but are largely absorbed in the wave function renormalization when passing to 
an effective theory of the low energy excitations. The issue here is whether there are 
correlations between (tentative) quasi particles, surviving at arbitrarily low energy scale. 
Of course the above refers to these latter correlations only. 

It was important to distinguish various scales of distances on the k-lattice of momenta, 
to obtain a clear picture of the rather singular correlations in k-space, and to relate 
Anderson's [IJ results for the phase shift to results from the T-matrix approximation. 
Anderson's finite phase shift, calculated for two particles residing on the very same point 
of the k-lattice, was seen to be not an artefact of this special choice, but represents the 
generic behavior in the limit k, k' — > dT with k' — > k. This behavior is in fact correctly 
signalled by the corresponding limit of the phase angle of the scattering vertex calculated 
in T-matrix approximation. 

Two fermions on the Fermi surface of a two-dimensional system repel each other at very 
short distances in momentum space. We note that this "partial exclusion principle" [1] 
is not only "partial", but also less persistent than genuine statistical correlations such as 
Pauli exclusion: The expectation value (nk F t n k F j.) will rise when more and more particles 
are added to the system, and may come arbitrarily close to one. The amplitude for two 
particles on the same Fermi point is only partitioned among different eigenstates with 
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energies in an infinitesimal interval at 2e^. However, the independence of quasi particles is 
obviously spoiled by these correlations at short distances in k-space. Since quasi particles 
on different Fermi points are however uncorrelated in 2D, Landau parameters involving a 
smooth angular average may still be well defined. 

The orthogonality effect caused by addition of an extra particle on the Fermi surface 
of a two-dimensional Fermi gas was shown to be finite within a crude approximation 
which takes into account only two-particle correlations between the extra particle and 
other particles in the system. This confirms earlier consistency checks of Fermi liquid 
theory performed directly in the infinite volume limit, J7|, || and a recent study showing 
the irrelevance of finite size effects in the T-matrix approximation. ]3[] In contrast to the 
case of addition of a local scatterer, which modifies the wave functions of other particles 
over the whole Fermi surface, a quasi particle in 2D modifies the wave functions at best 
on that Fermi point where it is added. In one dimension an orthogonality catastrophy 
does occur as a consequence of finite phase shifts for particles on opposite Fermi points. 

All well established weak coupling instabilities of the Fermi liquid are signalled by the 
renormalization group, evaluated perturbatively to some low order. |L4], [L7], [19|] A recent 
analysis of the crossover from ID Luttinger liquid behavior to 2D Fermi liquid behavior as 
a function of continuous dimensionality within perturbation theory summed to all orders 
indicates that at weak coupling higher orders in perturbation theory do not destroy the 
Fermi liquid fixed point in any dimension above one, as long as no Cooper instability 
sets in. [^0| In addition, recent rigorous results on two-dimensional Fermi systems seem to 
indicate that the existence of hitherto unknown weak coupling instabilities of the Fermi 



liquid is unlikely. [T^] These rigorous results are however not yet general enough to be 
applied to a system like the Hubbard model. 

On the other hand, the rather peculiar change of the two-particle wavefunction, espe- 
cially in 2D, could throw doubts on the general validity of conventional many-body theory 
itself, even if summed to all orders. Two particles near the same Fermi point indeed de- 
velop singular correlations in k-space, which are however visible only if the discrete fine 
structure of k-space is resolved. These correlations might not be adequately taken into 
account when taking the infinite volume limit before solving the full interacting prob- 
lem. In a two-dimensional Fermi liquid ground state, pairs of up- and down-spins near a 
common Fermi point would seem quite unstable objects: they have a tendency to repel 
each other but they can't, being blocked by their neighbors in k-space. The hypothesis 
of a complete non-perturbative reorganization of the ground state is therefore not com- 
pletely unplausible. In one dimension the repulsion of particles from common points in 
momentum space affects only excitations, not the ground state, which is instead modified 
by interactions between opposite Fermi points. However, in two dimensions the situa- 
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tion may be different, since new gapless degrees of freedom parallel to the Fermi surface 
appear. 

In summary, in our opinion there are interesting hints but no evidence for a breakdown 
of Fermi liquid theory at weak coupling in two dimensions. Clear is only that such a 
breakdown would have to be much more subtle that in one dimension. 
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Appendix A: Phase shift, phase angle and overlap in ID 

In one dimension, the non-interacting two-particle spectrum {E®} of states with fixed 
total momentum P is locally invariant, i.e. the levels become equidistant in the limit 
L — > oo, E® — > E, with a spacing AE® = E® +1 — E® that depends only on E. The level 
degeneracy d° a also becomes a function of E only. In this situation, the phase shift, energy 
shift and wave function overlap are uniquely determined by the phase angle, for almost 
all states in the large system limit, as we will now show. 

For large L, and assuming local spectral invariance, the eigenvalue equation (3.6) can 



be written as [H, 21 



jj — K l £ c (E) + ReK(E) (Al) 



where K(E) := lim^o lim^oo K L {E + irj), and 



d° 



is the "local sum", which is determined by levels in an infinitesimal interval centered at 
E. Here AE° = AE°(E) is the level spacing, and d° = d°(E) the degeneracy of levels 
with energy E. This decomposition holds asymptotically for almost all levels, exceptions 
being levels which are so close to the bottom or top of the spectrum that their local 
sum does not extend over a large number of levels on both sides. Using the identity 
J2m=-oo( x ~ m ) 1 = 7rcot(7rx), one obtains 

K l L ° c (E) = 7rp(E)cot(- X (E)) (A3) 

where p(E) = d°/L d AE° is the density of states, and 

X (E) := -irSE/AE (A4) 

the phase shift. Inserting (A. 3) into the eigenvalue equation, and solving for x(-E), one 
finds 

Unp(E) 



X (E) = - tan" 1 



(A5) 



1 - UReK(E) 

Since irp(E) = —ImK(E), this is nothing but the phase angle (f>(E), defined in (3.23). 
Hence, for locally invariant spectra, phase shift and phase angle are indeed equivalent for 
large systems. A slightly different derivation of this result can be found in Ref . f| . 
The normalization constant C is given by 
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for L — > oo and local spectral invariance. Note that here only levels in an infinitesimal 
neighborhood around E contribute. Using the identity J2m=-oo( x ~ m ) 2 — 71-2 '/ (sin(7r:r)) 2 , 
one obtains 

/ — I sin y 

C{E) = v^L^ (A.7) 
The overlap between interacting and non-interacting wave functions is thus 

(* a |tf a ) =: S(E) = Vd^C/5E = (A.8) 



Recall that is a symmetric combination of d° degenerate states with amplitude 1/vdP 
for each. The overlap of |$kk') and |^kk') is obtained from this by inserting E = + e^. 
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Figure Captions 



Fig. 1 Fermi surfaces dT and dT + on the discrete k-lattice in a finite system; the con- 
tinuous line represents the Fermi surface dT in the large system limit. 

Fig. 2 Regimes for the two-particle density of states in the (P, a;)-plane, separated by the 
functions u + (P), uj_(P) and ujq{P). The dotted line indicates values for (P,u) with 
p = x/4-71, x = 0.3. 

Fig. 3 Geometry of available two-particle states with fixed total momentum P and fixed 
energy in two dimensions. The bold sections on the circle around P/2 indicate 
the possible locations of momenta k and k' outside T such that k + k' = P and 
e k + e k' = const. 

Fig. 4 Two particles in k /m =k on the same Fermi point in 2D: generic case (a) and the two 
exceptional cases (b) and (c). Note that the plot shows an infinitesimal fraction of 
the Fermi surface, which therefore looks perfectly flat. 

Fig. 5 Phase shifts (points) and phase angle (line) as function of energy for a Hubbard 
model with next-neighbor hopping t — 1, coupling U — 5 and e F = —1.3 in ID, for 
a fixed total momentum P = 0. The phase shifts have been calculated for a finite 
system with L = 1000. 

Fig. 6 Phase shifts (dots), phase angle (solid line) and mean phase shift (dashed line) as 
function of energy for a Hubbard model with next- neighbor hopping t — 1, coupling 
U — 5 and e F = —2.3 in 2D, for a fixed total momentum P = (0.287T, 0.367r). The 
phase shifts have been calculated for a finite system with L = 100. Each point x{E) 
is obtained by averaging 50 phase shifts corresponding to eigenenergies near E (the 
dashed line connects these points as guide to the eye). 

Fig. 7 Phase angle for two particles with momenta k and k' close to a common point of 
the Fermi surface, plotted for various fixed total momenta P = k + k' as a function 
of the ratio |k' — k|/(P — 2kp ). A quadratic with a cutoff u c = 10, and a constant 
coupling U — 5 has been used in this plot. 
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